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The LiHoa;Yi_a;F4 Ising magnetic material subject to a magnetic field, Bx, perpendicular to the 
Ho''"'" Ising direction has shown over the past twenty years to be a host of very interesting ther- 
modynamic and magnetic phenomena. Unfortunately, the availability of other magnetic materials 
other than LiHoa;Yi_a;F4 that may be described by a transverse field Ising model remains very 
much limited. It is in this context that we use here mean-field theory to investigate the suitability 
of the Ho(OH)3, Dy(OH)3 and Tb(0H)3 insulating hexagonal dipolar Ising-like ferromagnets for 
the study of the quantum phase transition induced by a magnetic field, Bx, applied perpendicular 
to the Ising spin direction. Experimentally, the zero field critical (Curie) temperatures are known to 
be « 2.54 K, ^ 3.48 K and fa 3.72 K, for Ho(OH)3, Dy(0H)3 and Tb(OH)3, respectively. 
From our calculations we estimate the critical transverse field, B^, to destroy ferromagnetic order 
at zero temperature to be =4.35 T, Bl =5.03 T and Bl =54.81 T for Ho(OH)3, Dy(0H)3 and 
Tb(0H)3, respectively. We find that Ho(OH)3, similarly to LiHoF4, can be quantitatively described 
by an effective S = 1/2 transverse field Ising model (TFIM). This is not the case for Dy(OH)3 due 
to the strong admixing between the ground doublet and first excited doublet induced by the dipolar 
interactions. Furthermore, we find that the paramagnetic (PM) to ferromagnetic (FM) transition 
in Dy(OH)3 becomes first order for strong Bx and low temperatures. Hence, the PM to FM zero 
temperature transition in Dy(OH)3 may be first order and not quantum critical. We investigate 
the effect of competing antiferromagnetic nearest-neighbor exchange and applied magnetic field Bz 
along the Ising spin direction z on the first order transition in Dy(0H)3. We conclude from these pre- 
liminary calculations that Ho(OH)3 and Dy(OH)3, and their Y'^"^ diamagnetically diluted variants, 
Hoa;Yi_i(OH)3 and Dya;Yi_i(OH)3, are potentially interesting systems to study transverse-field 
induced quantum fluctuations effects in hard axis (Ising-like) magnetic materials. 



I. INTRODUCTION 
A. Transverse field Ising model 

Quantum phase transitions occur near zero tempera- 
ture and are driven by quantum mechanical fluctuations 
associated with the Heisenberg uncertainty principle and 
not by thermal fluctuations as in the case of classical 
temperature-driven phase transitions^*^. There is accu- 
mulating evidence that the exotic behavior exhibited by 
several metallic, magnetic and superconducting materials 
may have its origin in underlying large quantum fluctua- 
tions and proximity to a quantum phase transition. For 
this reason, much efforts are currently being devoted to 
understand quantum phase transitions in a wide variety 
of condensed matter systems. 

Perhaps the simplest model that embodies the phe- 
nomenon of a quantum phase transition is the transverse 
flcld Ising model (TFIM)— i^, first proposed by de Gennes 
to describe proton tunneling in ferrolectric materials^. 
The spin Hamiltonian for the TFIM reads: 

i^TFiM = E M\rj-r,)\)St{n)SI{r,) (1) 

(r-i) 



The Si {S ~ 1/2) quantum spin operators reside on the 
lattice sites r; of some d— dimensional lattice. The Si op- 
erators are related to the Pauli spin matrices by Si = ^CTi 
(here we set /i = 1). The components of Si obey the com- 
mutation relations [5f , S^] = i^ap-ySj 8ij where a, [3 and 
7 indicate x, y and z spin components, 5ij and eap-y are 
the Kronecker delta and fully antisymmetric tensor, re- 
spectively, r is the effective transverse field along the x 
direction, perpendicular to the Ising z axis. As described 
below, the effective transverse field T does not correspond 
one-to-one to the applied physical field Bx- Rather, T in 
Eq. ([2]) is a function of the real physical field B^- In 
fact, Jij is also a function of B^ If the spin inter- 

actions, Jij, possess translational invariance, the system 
displays for F = conventional long range magnetic order 
below some critical temperature, T,.. In the simplest sce- 
nario, where Jij > 0, the ordered phase is ferromagnetic 
and the order parameter is the average magnetization 
per spin, = {^/N)^^{Sf), where N is the number 
of spins. Since S'f and Sf do not commute, nonzero F 
causes quantum tunneling between the spin- up, |t), and 
spin-down, ||), eigenstates of Sf. By increasing F, Tc 
decreases until, ultimately, Tc{T) vanishes at a quantum 
critical point where F = Fc. On the T = tempera- 
ture axis, the system is in a long-range ordered phase 
for F < Fc while it is in a quantum paramagnetic phase 
for F > Fc. The phase transition between the paramag- 
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netic and long range ordered phase at constitutes the 
quantum phase transition^. 

One can also consider generalizations of the H^fim 
where the are quenched (frozen) random interactions. 
Of particular interest is the situation where there are 
as many ferromagnetic Jy > and antiferromagnetic 
Jij < couplings. This causes a high level of random 
frustration and the system, provided it is three dimen- 
sional, freezes into a spin glass state via a true thermo- 
dynamic phase transition at a spin glass critical temper- 
ature TjpiiiS. Here as well, one can investigate how the 
spin glass transition is affected by a transverse field F. 
As in the previous example, Tg(T) decreases as F is in- 
creased from zero until, at F = Fc, a quantum phase tran- 
sition between a quantum paramagnet and a spin glass 
phase ensues. Extensive numerical studies have found 
the quantum phase transition between a quantum para- 
magnet and a spin glass phas o^^i^^'^"^ to be quite interest- 
ing due to the occurrence of Griffiths-McCoy singularities 
(GMS)i^'i^. These GMS arise from rare spatial regions 
of disorder which may, for example, resemble the other- 
wise non-random (disorder-free) version of the system at 
stake. As a result, GMS can lead to singularities in var- 
ious thermodynamic quantities away from the quantum 
critical point. 



B. LiHo,Yi-^,F4 

On the experimental side, most studies 
aimed at exploring the phenomena associated 
with the TFIM have focused on the insulat- 
ing LiHo2;Yi_a;F4 (LHYF) Ising magnetic mate- 

j.^j^ ^l6. 17. 18. 19. 20. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32 this 

system, the Ho^+ Ising spin direction is parallel to the c 
axis of the body-centered tetragonal structure of LHYF. 
The random disorder is introduced by diluting the 
magnetic Ho'^"'" ions by non-magnetic Y'^"''. Crystal field 
effects lift the degeneracy of the ^Ig electronic manifold, 
giving an Ising ground doublet and a first excited 

|3>e) singlet at approximately 11 K above the ground 
doubleii^^. The other 14 crystal field states lie at much 
higher energiesS^. Quantum spin flip fluctuations are 
introduced by the application of a magnetic field, B^, 
perpendicular to the Ising c axis. admixes |$e) with 
|<I>j^), splitting the latter and producing an effective 
TFIM with r{B^)(x B^ for small B^ 

The properties of pure LiHoF4 in a transverse B^ are 
now generally qualitatively well understood^. Indeed, a 
recent quantum Monte Carlo study^ found general agree- 
ment between experiments and a microscopic model of 
LiHoF4. However, some quantitative discrepancies be- 
tween Monte Carlo and experimental data, even near the 
classical paramagnetic to ferromagnetic transition where 
Bx/Tc is small, do exist— i^. One noteworthy effect at play 
in LHYF at low temperatures is the significant enhance- 
ment of the zero temperature critical B^, B^, caused by 
the strong hyperfine nuclear interactions in Ho"^+-bascd 



materials^i^i^. 

LiHo2;Yi_a;F4 in a transverse B^ and a; < 1 has long 
been known to display paradoxical behaviors, both in 
the ferromagnetic (FM) (0.25 < a; < 1.0) and spin glass 
(SG) {x < 0.25) regimes. In the FM regime, a mean- 
field behavior Tc{x) oc x for the PM to FM transition 
is observed when B^ = 0—. However, in nonzero B^-, 
the rate at which Tc[Bx) is reduced hy B^ > in- 
creases faster than mean-field theory predicts as x is re- 
duced2^i21. In the high Ho'^'*' (SG) dilution regime (e.g. 
LiHoo.i67Yo.833F4), LHYF has long hcex^^l^^^^ argued 
to display a conventional SG transition for B^ = sig- 
nalled by a nonlinear magnetic susceptibility, X3 ^ diverg- 
ing at Tg as xz{T) cx (T - Tg)-^ iS. However, xz{T) 
becomes less singular as B^ is increased from B^. 
suggesting that no quantum phase transition between 
a PM and a SG state exists as T -i^"^. Recent 
theoretical studies^ i"^^'^^ suggest that for dipole-coupled 
Ho'^"'' in a diluted sample, nonzero B^ generates lon- 
gitudinal (along the Ising z direction) random fields 
that (i) lead to a faster decrease of Tc{Bx) in the FM 
regim022i21i2^ and (n) destroy the PM to SG transition 
for samples that otherwise show a SG transition when 

Bx = 0808.19.29^^^^^ oj.^ g^^. igj^g^^ ig^j to ^ Jigg,p_ 

pearance of the X3 divergence as B^ is increased from 



Perhaps most interesting among the phenomena ex- 
hibited by LHYF is the one referred to as anti- 
glass and which has been predominantly investigated in 
LiHoo.o45Yo.955F4^'23.24^^^ j.^^ reason for this name 

comes from AC susceptibility data on LiHoo.o45Yo.955F4 
which show that the distribution of relaxation times nar- 
rows upon cooling below 300 mK 16,23,24 ^ This behavior 
is quite different from that observed in conventional spin 
glasses where the distribution of relaxation times broad- 
ens upon approaching a spin glass transition at Tg > 
^il^. The antiglass behavior has been interpreted as evi- 
dence that the spin glass transition in LiHoa:Yi_a;F4 dis- 
appears at some nonzero Xc > 0. Results from more 
recent experimental studies on LiHoo.i65Yo.835F4 {x = 
16.5%) and LiHoo.o45Yo.955F4 {x = 4.5%) suggest an ab- 
sence of a genuine spin glass transition, even for a con- 
centration of Ho as large as 16.5%— In particular, it 
is in stark contrast with theoretical arguments^ which 
predict that, because of the long-ranged 1/r'^ nature of 
dipolar interactions, classical dipolar Ising spin glasses 
should have Tg{x) > for all a; > 0. However, even more 
recent work asserts that there is indeed a thermodynamic 
SG transition for x = 16.5% — , but that the behavior 
found in LiHoo.o45Yo.955F4 is truly unconventional^. 

Two very different scenarios for the failure of 
LiHoo.o4.5Yo.955F4 to show a spin glass transition have 
been put forwar d^^'^^i'^^'^° . Firstly, it has been suggested 
that the (small) off-diagonal part of the dipolar inter- 
actions lead to virtual crystal field excitations that ad- 
mix |tf>J) with |<I>e) and give rise to non-magnetic sin- 
glets for spatially close pairs of Ho'^"'" ions. The forma- 
tion of these singlets would thwart the development of 
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a spin glass state. This mechanism is analogous to the 
one leading to the formation of the random singlet state 
in dilute antiferromagnetically coupled 5 = 1/2 Heisen- 
berg spins^. However, a recent studj*^ shows that the 
energy scale for this singlet formation is very low (~ 10° 
mK) and that the random singlet mechanism^l may not 
be very effective at destroying the spin glass state in 
LiHoo.o45Yo.955F4 — . Hence the proposed formation of 
an entangled state in LiHoo.o45Yo.955F4 may, if it really 
exist, perhaps proceed via a more complex scheme than 
that proposed in Ref. Also, the low-temperature 

features observed in the specific heat in Ref. [2J| have 
not been observed in a more recent study^. Secondly, 
and from a completely different perspective, numerical 
simulations of classical Ising dipoles found that the spin 
glass transition temperature, Tg appears to vanish for a 
concentration of dipoles below approximately 20% of the 
sites occupied^i^i^. However, even more recent Monte 
Carlo simulations find that this conclusion may not be 
that firmly established^. 

As another possible and yet unexplored scenario, we 
note here that since Ho'^"'' is an even electron system 
(i.e. a non-Kramers ion), the Kramers' theorem is in- 
operative and the ground state doublet can be split 
by random (electrostatic) crystal field effects that com- 
pete with the collective spin glass behavior. For exam- 
ple, random strains, which may come from the substi- 
tution of Ho'^+— >Y'^+, break the local tetragonal sym- 
metry and introduces (random) crystal fields operators 
(e.g. Oj^) which have nonzero matrix elements be- 
tween the two states \^q) with |<I>g ) of the ground 
doublet, splitting it, and possibly destroying the spin 
glass phase at low Ho'^+ concentration. Indeed, such 
random transverse fields have been identified in sam- 
ples with very dilute Ho'^"'" in a LiYF4 matri x^^i^^ . 
Also, very weak random strains, hence effective ran- 
dom transverse fields, arise from the different (ran- 
dom) anharmonic zero point motion of ^Li and ^Li in 
Ho:LiYF4 samples with natural abundance of ^Li and 
^Li — . Finally, there may be intrinsic strains in the 
crystalline samples that do not arise from the Ho'^+/Y'^+ 
or ^Li/^Li admisture^. However, using available esti- 
mates^i^i^, calculations suggests that strain-induced 
random fields at play in LiHoo.o45Yo.955F4 may be too 
small [< 0(10^ mK)] to cause the destruction of the spin 
glass phase in this system^. Nevertheless, the point re- 
mains that, in principle, the non-Kramers nature of Ho'^"^ 
does offer a route for the destruction of the spin glass 
phase in LiHo2;Yi_a;F4 outside strictly pairwise, quan- 
tum^id^ or classical^2i2Sd£, magnetic interaction mecha- 
nisms. At this stage, this is clearly a matter that needs 
to be investigated experimentally further. One notes 
that, because of Kramers' theorem, the destruction of a 
SG phase via strain-induced effective random transverse 
fields would not occur for an odd-electron (Kramers) ion 
such as Dy'^+ or Er^+. In that context, one might think 
that a comparison of the behavior of LiDy3:Yi_a;F4 or 
LiEr2;Yi_j;F4 with that of LiHoi:Yi_2:F4 would be inter- 



esting. Unfortunately, while LiH0:cYi_2;F4 is an Ising 
system, the Er'^+ and Dy'^^ moments in LiEr2:Yi_2:F4 
and LiDyj.Yi„2:F4 are XY-like^i^. Hence, one can- 
not compare the LiEr3;Yi_3;F4 and LiDyi:Yi_2:F4 XY 
compounds with the LiHo2-Yi_2:F4 Ising material on the 
same footing. 

From the above discussion, it is clear that there are a 
number of fundamental questions raised by experimen- 
tal studies of LiHo2:Yi_2:F4, both in zero and nonzero 
transverse field B^, that warrant systematic experimen- 
tal investigations in other similar diamagnctically-diluted 
dipolar Ising-likc magnetic materials. Specific questions 
are: 

1. How does the quantum criticality of a transverse 
field Ising material with much smaller hyperfine in- 
teractions than Ho^^^ in LiHoa;Yi_i:F4 manifest it- 

2. Is the theoretical proposal of transverse- induced 
random longitudinal fields in diluted dipolar Ising 
materials^i^i^ valid and can it be explored 
and confirmed in other materials other than 
LiHo2:Yi_a;F^? In particular, are the phenom- 
ena observed in Ref. [27[ and ascribed to Griffiths 
singularities observed in other disordered dipolar 
Ising systems subject to a transverse field? 

3. Does the antiglass phcnomcnoni^i^iSi occur in 
other diluted dipolar Ising materials? If yes, which, 
if any, of the aforementioned theoretical proposals 
for the destruction of the spin glass state at small 
spin concentration is correct? 

C. RE(OH)3 materials 

As mentioned above, these questions cannot be inves- 
tigated with the LiEra;Yi_a;F4 and LiDya;Yi_2;F4 mate- 
rials isotructural to the LiHoa;Yi_a;F4 Ising compound 
since they are XY-like systems. However, we note in 
passing that it would nevertheless be interesting to ex- 
plore the topic of induced random fields^i^i^ and the 
possible existence of an XY dipolar spin glass and/or 
antiglass state in LiEr^Yi_^F4 and LiDya:Yi_a;F4. The 
LiTbj;Yi_a;F4 material is of limited use in such investi- 
gations since the single ion ground state of Tb'^+ in this 
compound consists of two separated singlets^, and lo- 
cal moment magnetism on the Tb'^^ site disappears at 
low Tb concentration^. In this paper, we propose that 
the RE(0H)3 (RE=Ho, Dy) compounds may offer them- 
selves as an attractive class of materials to study the 
above questions. Similarly to the LiHoF4, the RE(0H)3 
materials possess the following interesting properties: 

1. They are insulating rare earth materials. 

2. Their main spin-spin couplings are magnetostatic 
dipolc-dipole interactions. 
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3. The RE(0H)3 
perature. 



materials are stable at room tem- 



4. Both pure RE(0H)3 and LiHoF4 are colinear 
(Ising-like) dipolar ferromagnets with the Ising di- 
rection along the c axis of a hexagonal unit cell 
(RE(0H)3) or body-centered tetragonal unit cell 
(LiHoF4). In both cases there are two magnetically 
equivalent ions per unit cell. 

5. In RE(0H)3, the Kramers (Dy3+) and non- 
Kramers (Ho'^"'') variants possess a common crys- 
talline structure and both have similar bulk mag- 
netic properties in zero transverse magnetic field 
Bx- 

6. The critical temperature of the pure RE(0H)3 com- 
pounds is relatively high, ~ 3 K. This would make 
possible the study of Y— substituted Dy and Ho 
hydroxides down to quite low concentration of rare 
earth while maintaining the relevant magnetic tem- 
perature scale above the lowest attainable temper- 
ature with a commercial dilution refrigerator. 

7. Finally, and this is a key feature that motivated the 
present study, the first excited crystal field state 
in the Ho(OH)3 and Dy(0H)3 compounds is low- 
lying, hence allowing a possible transverse-field in- 
duced admixing and, possibly, a transverse field 
Ising model descriptior 
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To the best of our knowledge, it appears that the 
RE(0H)3 materials have so far not been investigated as 
potential realization of the TFIM. The purpose of this 
paper is to explore (i) the possible description of these 
materials as a TFIM, (ii) obtain an estimate of what the 
zero temperature critical transverse field may be and, 
(iii) assess if any new interesting phenomenology may oc- 
cur, even in the pure compounds, in nonzero transverse 
field Bx. 

We note, however, that there are so far no very large 
single crystals of RE(0II)3 available^. For example, 
their length typically varies between 3 mm and 17 mm 
and their diameter between 0.2 and 0.6 mm. The lack of 
large single crystals would make difficult neutron scat- 
tering experiments. However, possibly motivated by 
this work and by a first generation of bulk measure- 
ments (e.g. susceptibility, specific heat), experimental- 
ists and solid state chemists may be able to conceive 
ways to grow larger single crystals of RE(0H)3. Also, in 
light of the fact that most experiments on LiHo2.Yi_2:F4 
that have revealed exotic behavior are bulk measure- 
mentsi^iiiii^i^iS^i^i^, we hope that at this time the lack 
of availability of large single crystals of the RE(0H)3 se- 
ries is not a strong impediment against pursuing a first 
generation of bulk experiments on RE(0H)3. 

The rest of this paper is organized as follows. In Sec- 
tion II, we review the main, single ion, magnetic prop- 
erties of RE(0H)3 (RE=Dy, Ho, Tb). In particular, we 
discuss the crystal field Hamiltonian of these materials 







O'-, y 




y 


Tb(0H)3 


0.3952(7) 


0.3120(6) 


0.276(1) 


0.142(1) 


Dy(OH)3 


0.3947(6) 


0.3109(6) 


0.29(3) 


0.15(2) 


Ho(OH)3 


0.3951(7) 


0.3112(7) 


0.30(3) 


0.17(3) 


Y(0H)3 


0.3958(6) 


0.3116(6) 


0.28(1) 


0.17(1) 



Table I: Position parameters of O ~ and H~ ions in rare earth 
hydroxides and Y(0H)3 from Ref. [5^ and Ref. [5^ (see text 
in Section Uraj. 



and the dependence of the low-lying crystal field levels 
on an applied transverse field Bx- We present in Sec- 
tion III a mean-field calculation to estimate the Bx vs 
temperature, T, Tc{Bx) phase diagram of these materi- 
als. In Section IV, we show that Ho(OH)3 and Tb(0H)3 
can be described quantitatively well by a transverse field 
Ising model, while Dy(0H)3 cannot. The Subsection A 
of Section V, uses a Ginzburg-Landau theory to explore 
the first order paramagnetic (PM) to ferromagnetic (FM) 
transition that occurs in Dy(0H)3 at low temperatures 
and strong Bx- The following Subsection B discusses 
the effect of nearest-neighbor antiferromagnetic exchange 
interaction and applied longitudinal (i.e. along the z 
axis) magnetic field, Bz, on the first order transition in 
Dy(0H)3. A brief conclusion is presented in Section VI. 
Appendix A discusses how the excited crystal field states 
in Dy(0H)3 play an important quantitative role on the 
determination of Tc{Bx) in this material. 



II. RE(OH)3: MATERIAL PROPERTIES 



A. Crystal properties 



All the rare earth hydroxides form hexagonal crys- 
tals that that are iso-structural with Y(0H)3. The lat- 
tice is described by translation vectors ai = (0,0,0), 
a2 = {—a/2, 0) and as (0, 0, c). A unit cell con- 

sist of two Ho'^+ ions at coordinates (1/3,2/3,1/4) and 
(2/3, 1/3,3/4) in the basis of lattice vectors Oi, 02 and 
as- The coordinates of three O^" and H^ ions, relative 
to the position of Ho-^+j are ±{x,y,0), ±{—y,x — y,0) 
and — X, —x,0), where ± refers to the first and sec- 
ond Ho'^+ in the unit cell, respectively^. The values of 
the parameters x and y are listed in Table IH x=0.396, 
y=0.312 and for H": x=0.28, y=0.17^. The lattice 
structure is depicted in Fig. [TJ The lattice constants for 
Tb(0H)3, Dy(0H)3, Ho(OH)3 and Y(0H)3, from Beall 
et al<^, are collected in Table [Hi Each magnetic ion is 
surrounded by 9 oxygen atoms that create a crystalline 
field characterized by the point group symmetry C^ir^. 
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a 


c 


c/a 


Tb(0H)3 


6.315(4) 


3.603(2) 


0.570(5) 


Dy(0H)3 


6.286(3) 


3.577(1) 


0.569(0) 


Ho(OH)3 


6.266(2) 


3.553(1) 


0.567(0) 


Y(0H)3 


6.261(2) 


3.544(1) 


0.566(0) 



Table II: Lattice constants for rare earth hydroxides and 
Y(OH)3 (from Ref. 




Figure 1: The lattice structure of rare earth hydroxides and 
Y(OH)3. The arrows indicate the 6 nearest neighbors of the 
central ion in the lower prism. The lower left corner shows 
a projection of the unit cell onto the plane given by lattice 
vectors ai and a^- 



B. Single ion properties 

The electronic configuration of the magnetic ions 
Tb3+, Dy3+ and }io^+ is, respectively, 4/^, 4/^ and 4/i". 
Magnetic properties of the rare earth ions can be de- 
scribed by the states of the lowest energy multiplet: the 
spin-obit splitting between the ground state J manifold 
and the first excited states is of order of few thousands 
K. The ground state manifolds can be found from Hund's 
rules and are '^Fq, and ^/s for Tb'''+, Dy^+ and 

Ho'^"'", respectively. The Wigner-Eckart theorem gives 
the Lande (^-factor equal to |, | and | for Tb'^+, Dy'^+ 
and Ho^+, respectively. 

In a crystalline environment, an ion is subject to the 
electric field and covalency effects from the surrounding 
ions. This crystalline field effect partially lifts the de- 
generacy of the ground state multiplet. The low energy 
levels of Tb'^+ in Tb(0H)3 are a pair of singlets, that 
consist of the symmetric combination of |±6) states with 
a small admixture of the |0) state and an antisymmetric 
combination 0.3 cm~^ above^. The next excited state is 
well separated from the lowest energy pair by an energy 
of 118 cm-i ^ (1 cm-i « 1.44 K). In the case of Dy3+ 
in Dy(0H)3, the spectrum consist of 8 Kramers doublets 
with the first excited state 7.8 cm~^ above the ground 



is composed of a ground state doublet and an excited 
singlet state 11.1 cm"^ above^. 

Due to the strong shielding of the 4/ electrons by the 
electrons of the filled outer electronic shells, the exchange 
interactions for 4/ electrons is weak and the crystal field 
can be considered as a perturbation to the fixed J man- 
ifold. Furthermore, because the strong spin-obit interac- 
tion yields a large energy gap between the ground state 
multiplet and the excited levels, we neglect all the excited 
electronic multiplets in the calculation. 

According to arguments provided by Stevens^, we ex- 
press the matrix elements of the crystal field Hamiltonian 
for the ground state manifold in terms of operator equiv- 
alents. The details of the method and conventions for 
expressing the crystal field Hamiltonian can be found in 
the review by Hutchings^. On the basis of the Wigner- 
Eckart theorem, one can write the crystal field Hamilto- 
nian in the form 



(2) 



where O™ are Steven's "operator equivalents". On are 
constants called Stevens multiplicative factors and B"!^ 
are crystal field parameters (CFP). The CFP are usually 
determined by fitting experimental (spectroscopic) data. 
From angular momentum algebra, we know that in the 
case of f electrons, we need to consider only n = 0, 2, 4, 6 
in the sum ([2]). The choice of coefficients in Hamil- 
tonian ^ that do not vanish and have nonzero corre- 
sponding matrix elements is dictated by the point sym- 
metry group of the crystalline environment. The Stevens 
operators, O™, are conveniently expressed in terms of 
vector components of angular momentum operator J. In 
the case of the RE(0H)3 materials, considered herein, 
the point-symmetry group is C^h, and the crystal field 
Hamiltonian is of the form 



a.,BlOl{J,) + PjBlOl{J,) 
+ljBlOl{J,) + ^jBlOliJ,). (3) 



stated. The low energy spectrum of Ho^+ in Ho(OH)3 



The Stevens multiplicative factors a,/, (3,j and 7j (02, 6i 
and 6*6) are collected in Table. IIIII 

For the sake of conciseness, and to illustrate the proce- 
dure, most of our numerical results below are presented 
for one set of CFP only. The qualitative picture that 
emerges from our calculatiosn does not depend on the 
specific choice of CFP parameters. Only quantitative 
differences arc found using the different sets of CFP. Ulti- 
mately, a further experimental determination of accurate 
-B™ values would need to be carried out in order to obtain 
more precise mean-field estimates as well as to perform 
quantum Monte Carlo simulations of the Re(0H)3 sys- 
tems. According to our arbitrary choice^, if not stated 
otherwise, we use in the calculations the CFP provided 
by Scott et alM-^ll^^. For Ho(OH)3 and Dy(0H)3 
different values of CFP were proposed by Karmakar et 
al^^. As one can see in Fig. [31 for Ho(OH)3, the latter 
set of CFP yields a somewhat higher mean-field critical 
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Ion 


oj 


Pj 


7J 


Tb3+ 


-1/(32-11) 


2/(33 -5-112) 


-1/(34 ■ 7 - ll2 ■ 13) 


Dy3+ 


-2/(32 -5-7) 


-23/(33 -5-7-11-13) 


22/(33-7-112-132) 


Ho3+ 


-1/(2-32.52) 


-1/(2 - 3 - 5 - 7 - 11 - 13) 


-5/(33 . 7 . ii2 . 132) 



Table III: Stevens multiplicative factors 



Ref 


Crystal 


B2"(cm-i) 


B2(cm-') 


Bg(cm-i) 


B|(cm-i) 


54 


Tb(OH)3* 


207.9 ± 2.8 


-69.0 ± 1.6 


-45.3 ± 1.1 


585 ± 10 


54 


Tb:Y(OH)3 


189.1 ± 2.6 


-69.1 ± 1.5 


-45.7 ± 1.0 


606 ±9 


54 


Dy(OH)3* 


209.4 ± 3.4 


-75.5 ±3.5 


-40.1 ± 1.9 


541.8 ±26.5 


61 


Ho(OH)3 


200 ± 2.0 


-57 ±0.5 


-40 ±0.5 


400 ±5 


54 


Ho:Y(OH)3* 


246.0 ± 3.4 


-56.7 ± 1.2 


-39.8 ±0.3 


543.6 ± 3.3 


62 


Dy(OH)3 


215.9 


-72.0 


-42.0 


515.3 



Table IV: Crystal field parameters. Some of the calculation we performed with one set of crystal field parameters only. Crystal 
field parameters arbitrary chosen^ to be used in these calculations are marked with *. 



temperature and quite a bit higher critical value of the 
transverse magnetic field B'^ =7.35 T, compared with 
B% =4.35 T obtained using Scott et aVs CFP^(see Fig. 
E]). Similarly, Karmakar et aUs CFP-^ for Dy(0H)3 give 
a much higher critical field of =9.12 T, compared with 
= 5.03 T when Scott 's et al.'s CFF^i^i^ are used. 
From the two sets of CFP for Tb(0H)3 we choose the one 
obtained from measurements on pure Tb(0H)3^. Using 
the CFP obtained for the system with a dilute concen- 
tration of Tb in a Y(0H)3 matrix, Tb:Y(0H)3^, makes 
only a small change to the value of critical transverse 
field; we obtained B^ =50.0 T and B^ =54.8 T calcu- 
lated using Tb:Y(0H)3 and Tb(0H)3 CFP, respectively 
(see Fig. [3)). Available values of the CFP are given in 
Table El 

We show in Table fVl the lowest eigenstates and eigen- 
values of the crystal field Hamiltonian ([3]). The calcu- 
lated energies are not in full agreement with the experi- 
mentally determined values because the CFP were fitted 
using all the observed optical transitions, including tran- 
sitions between diffent J manifolds^. Furthermore, the 
fitting procedure used by Scott^ includes perturbative 
admixing between manifolds with the admixing incorpo- 
rated into effective Stevens multiplicative factors aj, pj 
and 7j that slightly differ from those given in Table IIIII 

Given the uncertainty in the CFP, which ultimately 
lead to an uncertainty of approximately ~ 40% on B^ 
in Ho(OH)3 and Dy(0H)3, as well as the nature of the 
mean-field calculations that we use and which neglects 
thermal and quantum mechanical fluctuations, as well as 
for simplicity sake, we ignore here the effect of hyperfine 
coupling of the electronic and nuclear magnetic moments. 
However, as shown for LiHoa;Yi_;;;F4, the important role 
of hyperfine interactions for Ho^+ on the precise deter- 
mination of B^ must eventually be considere d^'^^i^^'^"^ . 
At this time, one must await results from further exper- 
iments and a precise set of CFP for He in order to go 
beyond the mean-field calculations presented below or 



Eigenstate Energy [cm ^] 



Dy(OH)3 


0.92 |±15/2) - 0.15 |±3/2) + 0.37 |t9/2> 




0.40 |±15/2) ± 0.20 |±3/2) - 0.90 |=f9/2) 


9.6 


Dy(OH)3 Karmakar et al, Ref. [61^1 


0.98 1±15/2) - 0.09 |±3/2) + 0.15 |t9/2> 




0.17 I±15/2) ± 0.22 |±3/2) - 0.96 |t9/2> 


19.3 


Ho:Y(OH)3 


0.94 |±7) ± 0.31 1±1) ± 0.15 |t5) 




0.59 |6) ±0.55 |0) ±0.59 1-6) 


12.7 


Ho(OH)3 Karmakar et al., Ref. [61^1 


0.97 |±7) ± 0.24 |±1) ± 0.09 |=f5) 




0.60 |6) ±0.52 |0> ±0.60 |-6> 


23.6 


Tb(0H)3 


0.71 |6)± 0.05 |0)± 0.71 1-6) 




0.71 |6) -0.71|-6> 


0.49 


0.99 |±5) ±0.13 |=Fl) 


122.06 


Tb:Y(0H)3 


0.71 |6)± 0.06 |0>± 0.71 1-6) 




0.71 |6) -0.711-6) 


0.58 


0.99 |±5) ±0.16 |=Fl) 


115.33 



Table V: Eigenstates and energy levels calculated with the 
crystal field parameters collected in Table ITVl 



to pursue quantum Monte Carlo calculations as done in 
Refs. 00]. As suggested in Ref. [1], the accuracy of any 
future calculations (mean-field or quantum Monte Carlo) 
could be improved by the use of directly measured accu- 
rate values of the transverse field splitting of the ground 
state doublet instead of the less certain values calculated 
from CFP. 

Since our main goal in this exploratory work is to es- 
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timatc the critical transverse field, B^, for the family of 
RE(0H)3 compounds and to explore the possible vahd- 
ity of a transverse field Ising model description of these 
materials, we henceforth restrict ourselves to the Hd in 
Eq. ([2]) with the CFP (B™ parameter values) given in Ta- 
ble IIVI These calculations could be revisited and quan- 
tum Monte Carlo simulations^ii performed once experi- 
mental results reporting on the effect of on Dy(0H)3 
and Ho(OH)3 become available. 



C. Single ion transverse field spectrum 

A magnetic field, B^., applied in the direction trans- 
verse to the easy axis splits the degeneracy of the ground 
state doublet in the case of Ho(OH)3 and Dy(0H)3, or 
increase the separation of the ground levels in the case 
of the already weakly separated singlets in Tb(0H)3. By 
diagonalizing the single-ion Hamiltonian, Tio, which con- 
sist of the crystal field and Zeeman term, 



Tio = TicfiJi) - gi-J-BB^Jx 



(4) 



we obtain the transverse field dependence of the single ion 
energy levels, plotted in Fig. ^ In the case of Dy(0H)3, 
the two lowest energy levels splitting is too small to be 
clearly visible in the main panel of Fig. [21 Hence, we 
show the energy separation between the two lowest levels 
in the inset of Fig. [2] for Dy(0H)3. Furthermore, the 
separation vanishes at B^ = 3.92 T, indicating that the 
two lowest states for this specific value of the transverse 
field, Bx, are degenerate. 

To calculate the transverse field dependence of the low- 
est energy levels up to the critical transverse field where 
dipolar ferromagnetism is destroyed, we do not have to 
include all the crystal field states since the i?a;-induced 
admixing among the states decreases with increasing en- 
ergy separation. In the case of IIo(OII)3 we can repro- 
duce the field dependence, B^, of the lowest energy levels, 
E in Fig. [2] using only the four lowest levels. However, in 
order to achieve a similar level of agreement for Dy(0H)3, 
we have to retain the ground doublet and several of the 
lowest excited doublets. 



III. NUMERICAL SOLUTION 

The collective magnetic properties of the considered 
rare earth hydroxides are mainly controlled by a long 
range dipolar interaction between the magnetic moments 
carried by the rare earth ions. The dipolar interaction 
is complemented by a short range exchange interaction. 
Adding the interaction terms to the single ion Hamilto- 
nian (HI gives a full Hamiltonian, H, of the form 

Ti- = y^7Ycf(Ji) - gfJ-By^^BxJj^x 

i i 

+ 2 Ji,fj.jj-iy + 2'-^°^ 5Z ' 

ij i,nn 




30 



20 



10 



-15^ 



1 2 3 4 B [T] 



2 4 



-10 



-20 




Ho(OH)^ 



2 4 
B [T] 



Figure 2: Energy levels as a function of applied transverse 
field for Dy(OH)3 and Ho(OH)3 . The inset shows the sepa- 
ration of the lowest energy levels in Dy(OH)3. 




Figure 3: The main panel shows the phase diagrams for 
Dy(0H)3 and Ho:Y(OH)3 (crystal field parameters of Scott et 
^^_54j58'j rpj^g indicates the location of the tricritical point 
for Dy(OH)3. The transition is first order for temperatures 
below the temperature location of the tricritical point. The 
upper inset shows the phase diagram for Tb(OH)3 (crystal 



field parameters of Scott et aZ.- 



-). The lower inset shows 



the phase diagram for Ho(OH)3 and Dy(OH)3 calculated with 
the crystal field parameters of Karmakar et alM^^. 



C^!^ are the anisotropic dipole-dipolc interaction con- 
stants of the form C'^^ ~ ^"i^a'^'* ^ij^ where fijV ^ 
X, y, z; a is a lattice constant (see Tabic HI)) and fiQ is the 
permeability of vacuum. L'^^ are dimensionless dipolar 
interaction coefficients. 



(6) 
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where Tij = rj — Vi, with the lattice position of 
magnetic moment J^, expressed in units of the lattice 
constant, a. jTcx is the antiferromagnctic (jTcx > 0) 
exchange interaction constant, which can be recast as 

Je^ = ^°4ja^'' -^ex, whcrc Jex IS now a dimensionless 
exchange constant that, when multiplied by the nearest 
neighbor coordination number, z=6, can be use to com- 
pare the relative strength of exchange vs the magnetic 
dipolar lattice sum (energies) collected in Table IVTl The 
label nn in Eq. (O denotes the nearest neighbor sites of 
site i. 

The exchange interaction is expected to be of some- 
what lower strength than the dipolar coupling^i^. We 
therefore neglect it in most of the calculations, but we 
discuss its effect on the calculated vs T phase dia- 
gram at the end of this section as well as explore its in- 
fluence on the occurrence of a first order phase transition 
in Dy(0H)3 in Section IVB] Denoting L^^ = L'^^ and 

j^fj.v = ^"i^a^'* Lfj^i^, we write a mean-field Hamiltonian 
in the form 

Hmf ~ Ti-dJ) ~ gi-^BBxJx 

+ E (^MM + 2^ex) (jp (Jp> - ^ ( Jm)') (7) 
^—x,y.z ^ ^ 

with 2 = 6 the number of nearest neighbors. The last 
term in Eq. dH), — ^ {J^i)" , has no effect on the calculated 
thermal expectation values of the x and z components of 
the magnetization, and can be dropped. The off-diagonal 
terms, £p,y with =/= vanish due to the lattice symme- 
try. We employ the Ewald tcchnique^i^i^i^iia to cal- 
culate the dipole-dipole interaction, L^^ , of Eq. By 
summing over all sites j coupled to an arbitrary site i, 
we obtain the coefficients L^^ listed in Table IVII The 
considered Ewald sums ignore a demagnetization term^S 
and our calculations can therefore be interpreted as cor- 
responding to a long needle-shape sample. 





Lxx 






Tb(OH)3 


-11.43 


-11.43 


-28.01 


Dy(OH)3 


-11.40 


-11.41 


-28.20 


Ho(OH)3 


-11.38 


-11.37 


-28.45 



Table VI: Dimensionless lattice sums calculated with the val- 
ues oi c/a taken from Table HIl 

We diagonalize numerically TYmf in Eq. ([7]), and cal- 
culate self-consistently the thermal averages of Jx and Jz 
operators, from the expression 

Trace[J^ exp(-7YMF/?')] 
^ " Trace[exp(-7^MF/T)] ' ^ ' 

where ^ stands for x and z. (Jy) = due to the lattice 
mirror symmetries and since B is applied along x. 

For a given Bx, we find the value of the critical tem- 
perature, Tc{Bx), at which the order parameter, (Jz), 



vanishes. The resulting Bx vs T phase diagrams, ob- 
tained that way, using all sets of CFP from Table ITVl are 
shown in Fig. [3l In the main panel, we plot the phase 
diagrams for Ho(OH)3 and Dy(0H)3, using Scott et aZ.'s 
CFB^iSi58i64. The top inset shows the Bx vs T phase 
diagrams for Tb(0H)3, using two available sets of CFP. 
This indicates that, for Tb(0H)3, the critical field Bx{T) 
reaches very quickly the upper limit of magnetic fields 
attainable with commercial magnets. The bottom inset 
shows the Bx vs T phase phase diagrams for Ho(OH)3 
and Dy(0H)3 using Karmakar et aUs CFB^iS. Al- 
though the diagrams differ quantitatively for the two sets 
of CFP, the overall qualitative trend is the same for both 
sets. Table IVIII lists the mean-field estimates of Tc and 
B^ together with the experimental values of 7^ ^53,65,71 ^ 

There are two contributing factors behind the differ- 
ence between the experimental and mean-field values of 
Tc in Table IVIII and, presumably once they are experi- 
mentally determined, those for B^. Firstly, in obtaining 
those mean- field values from Eq. ([7|) and Eq. ([8]), we 
neglected the (presumably) antiferromagnteic nearest- 
neighbor exchange jTcx which would contribute to a de- 
pression of both the critical ferromagnetic temperature 
Tc and B^. Secondly, mean- field thgeory neglects corre- 
lations in the thermal and quantum fluctuations which 
would also contribute to reduce Tc and Bx- From the 
comparison of mean-field theory^ and quantum Monte 
Carlo^ for LiHoF4, we would anticipate that our mean- 
field estimates of Tc and Bx are accurate within 20% to 
40%, notwithstanding the uncertainty on the crystal field 
parameters. 

By seeking a self-consistent solution for {Jz), starting 
from either the fully polarized or weakly polarized state, 
two branches of solutions are obtained at low temper- 
ature and large Bx for Dy (011)3. This suggests a first 
order PM to FM transition when using cither set of CFP 
for this material. This result was confirmed by a more 
thorough investigation (see Section V below). The top 
right inset of Fig. [7] shows the behavior of the {Jz) as a 
function of Bx for T = 0.3 K, illustrating the transition 
field and the limits for the superheating and supercool- 
ing regime. The black dot in the main panel and inset 
of Fig. [3] shows the location of the tricritical point (see 
Section V) . Note that the Bx value at the tricritical point 
is ^ 4.85 T using the CFP of Scott et al. (main panel of 
Fig. 54, 57, 58, 64 ^ Hencc, the occurrence of a first order 
transition here is not directly connected to the degen- 
eracy occurring between the two lowest energy levels at 
Bx = 3.92 T using the same set of CFP (see inset of Fig. 
[5] for Dy(0II)3). A zoom on the low temperature regime 
and the vicinity of the tricritical point for Dy(0II)3 is 
shown in Fig. [71 The calculation details needed to ob- 
tain the phase diagram of Fig. [7] are described in Section 
IV Al The existence of a first order transition at strong Bx 
in Dy(0II)3 depends on the details of the chosen Hamil- 
tonian in Eq. For example, as discussed in Section 
IV Bl a sufficiently strong nearest-neighbor antiferromag- 
nctic exchange, jTcx; eliminates the first order transition. 
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Crystal | experimental Tc [K]|MFT Tc [K]|MFT Bl [T] 
Ho(OH)3 | 2.54 I 4.28 I 4.35 

Dy(0H)3 3^48 5£1 5.03 

Tb(OH)3| 3.72 I 5.59 | 54.81 

Table VII; Experimental values of critical temperatures 
and mean-field theory (MFT) estimates for Tc and 

B%. 

We also discuss in Section IV Bl the role of a longitudinal 
field Bz (along c axis) on the first order transition. At 
this time, one must await experimental results to ascer- 
tain the specific low temperature behavior that is at play 
for strong in Dy(0H)3. 

We now briefly analyze the effect of a nonzero exchange 
interaction. The dependence of the critical temperature, 
Tc, and the critical transverse field, -B^, on the exchange 
constant, Jex, is plotted in Fig. 21 The dot on the B^ vs 
Jex plot for Dy(0H)3 indicates the threshold value of Jex, 
J^^'^=0.995, above which the first order transition ceases 
to exist. The dependence of the existence of the first 
order transition on Jex is discussed in some detail in Sec- 
tion |VB] For Jex < Jex'^, the thinner lines correspond 
to the boundary of the supercooling and superheating 
regime. In the mean-field theory presented here, Jex sim- 
ply adds to the interaction constant L^^ with fi = x,y,z 
in Eq. ([7]) (see Table IVl)) . Hence, beyond a threshold 
value of Jex, the system no longer admits a long range 
ordered ferromagnetic phase. In the case of Dy (011)3, 
stays almost unchanged as a function of Jex, until it 
drops very sharply when Czz + 2i/ox = (z = 6). In the 
inset of Fig. [U we focus on the regime where B^ vs Jex 
plot sharply drops. The cusp at B^ =3.92 T is a conse- 
quence of the degeneracy of the lowest energy cigenstates 
(see Fig. As will be shown in detail in the next sec- 
tion for the Ho(OH)3 system, the energy gap separating 
transverse-field-splitted levels of the ground state dou- 
blet plays the role of an effective transverse field T{Bx) 
acting on effective Ising spins. 

IV. EFFECTIVE 5=1/2 HAMILTONIAN 

In this section we show that IIo(0H)3 and Tb(0II)3 
can be described with good accuracy by an effective 
TFIM Hamiltonian. On the other hand, although 
Dy(0H)3 has been referred in the literature as an Ising 
materia l^^i^^ , we find that it is not possible to describe 
the magnetic properties of this material within the frame- 
work of an effective Ising Hamiltonian that neglects the 
effect of the excited crystal field states. 

To be able to identify a material as a realization of 
an effective microscopic Ising model, the following con- 
ditions should applyi^: 

• There has to be a ground state doublet or a close 
pair of singlets that are separated from the next 




-1 012345 -1 012345 



Figure 4: The effect of the value of the exchange constant on 
the phase boundary, Tc vs Jex at Bx=Q (left) and B% vs Jox at 
T = (right). Solid and dashed lines are used for Dy(OH)3 
and Ho(OH)3, respectively. The inset of the right graph shows 
a focus on the features of the high Jex regime in the Dy(OH)3 
plot. The dot in the right panel indicates the tricritical point. 
The two additional lines at the left side of the tricritical point 
mark the limits of superheating and supercooling regimes. In 
the calculations, the CFP of Ref. [Hj were used. 

energy level by an energy gap that is large in com- 
parison with the critical temperature. This ensures 
that at the temperatures of interest only the two 
lowest levels are significantly populated. 

• To first order, there has to be no transverse sus- 
ceptibility. It means that there should be no ma- 
trix elements of (J^,, Jy) operators between the two 
states of the ground doublet. 

• Furthermore, the longitudinal (in the easy axis di- 
rection) susceptibility has to be predominantly con- 
trolled by the two lowest levels. In other words, 
there has to be no significant mixing of the states 
of the lowest doublet with the higher levels via the 
internal mean field along the Ising direction. In 
more technical terms, the van Vleek susceptibility 
should play a negligible role to the non-interacting 
(free ion) susceptibility near the critical tempera- 
ture2i. 

• In setting up the above conditions, one is in ef- 
fect requesting that a material be describable as a 
TFIM from a miscroscopic point of view. However, 
one can, alternatively, ask whether the quantum 
critical point of a given material is in the same uni- 
versality class as the relevant transverse field Ising 
model. In such long as transition is sec- 
ond order, then sufficiently close to the quantum 
critical point, a mapping to an effective TFIM is 
always in principle possible. However, it can be 
difficult to estimate the pertinent parameters enter- 
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ing the Ginzburg-Landau- Wilson theory describing 
the transition. From that perspective, the first and 
the third point above are always fulfilled sufficiently 
close to a second order quantum critical poin1*2^. 

The first condition is not satisfied in the case of 
Dy(0H)3. The energy gap of 7.8 cm'^ w 11.2 K is 
not much larger than the mean-field critical temperature 
Tc ^ 5.31 K. Hence, at temperatures close to Tc, the 
first excited doublet state is also significantly populated. 
Furthermore, and most importantly, in the context of a 
field-induced quantum phase transition, the third con- 
dition above is also not satisfied. Hence, even at low 
temperatures, because of the admixing of the two lowest 
energy states with the higher energy levels that is induced 
via the internal (mean) field from the surrounding ions, 
Dy(0H)3 cannot by described by an effective microscopic 
Ising model that solely considers the ground doublet and 
ignores the excited crystal field states. This effect and 
the associated role of nonzero Jz matrix elements be- 
tween the ground state and higher crystal field levels is 
discussed in more detail in Appendix [Xj As an interest- 
ing consequence of this participation of the higher energy 
levels, we predict that, unlike in the TFIM of Eq. (1), a 
first order phase transition may occur at high transverse 
field in Dy(0H)3 (see Section |V^. 

For Ho(OH)3 and Tb(0H)3 we construct an effective 
Ising Hamiltonian, following the method of Refs. [IQl!]. 
We diagonalize exactly the noninteracting Hamiltonian, 
Ho of Eq. (ID), for each value of the transverse field, Bx- 
We denote the two lowest states by \a{Bx)) and \(3{Bx)) 
and their energies by Ea{Bx) and Ei3{Bx), respectively, a 
transverse field enforces a unique choice of basis, in which 
the states can be interpreted as |— and |<— ) in the Ising 
subspace. We introduce a new |t) and ||) basis, in which 
the Jz matrix elements arc diagonal, by performing a 
rotation 

|T) = ^(|«(B,))+exp(i0)|/?(i?,))), 

U) = -l5(|<.(S.))-cxp(^0)|/3(i?.))). ^ ^ 

In this basis, the effective single ion Hamiltonian, de- 
scribing the two lowest states, is of the form 

Ht = EcuiBx) - ^A(B,)a^ (10) 

where EcuiB,) = \{E^{Bx) + Ep{Bx)) and A(S,) = 
Ep{Bx) — Ea{Bx)- Thus the splitting of the ground 



state doublet plays the role of a transverse magnetic field, 
r = iA(S^) in Eq. (1). In the case of Tb(0H)3, after 
performing the rotation even at B^ = 0, a small 
transverse field term (F = ^A(O) > 0) is present in 
Hamiltonian fTO]). For Dy(0H)3 and Ho(OH)3, the split- 
ting of the energy levels, obtained via exact diagonaliza- 
tion was already discussed at the end of Section [TTl and is 
shown in Fig. [21 To include the interaction terms in our 
Ising Hamiltonian, we expand the matrix elements of Jx , 
Jy and Jz operators in terms of the a'^ {v = x, y, z) Pauli 
matrices and a unit matrix, a'' = 1, 

J^^^^C,,1+ ^ C^,{BxK. (11) 
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Figure 5: Tlie coefficients in the projection of J operators 
onto tlie two dimensional Ising subspace for Ho(OH)3. The 
inset shows the ratios of the coefficients present in Eq. (|12p . 

By replacing all Ji^^ operators in the interaction term 
of Hamiltonian ^ by the two dimensional representation 
of Eq. pT|) . one obtains in general a lengthy Hamilto- 
nian containing all possible combinations of spin-i inter- 
actions. In the present case, the resulting Hamiltonian 
is considerably simplified by the crystal symmetries and 
the consequential vanishing of off-diagonal elements of 
the interaction matrix This would not be the case 

for diluted Hoa;Yi_2;F4 (see Ref. [1]]). After perform- 
ing the transformation in Eq. ([9|), we have J^ = Czzcf, 
Jf = Cyyal and Jf = C^x^^i + Cx^- Hence, we can 
rewrite the mean-field Hamiltonian ([7]) in the form 



Umf = {Czz + zJ^^)Cizmz(j'- + [CxxCxCxx ~ ^A(S^) ) 
+ [Cxx + zJc^)ClxmxCr'' + {Cyy 4- zJcyc)Cyymyay , 



(12) 



where rrii, = {a'^) and (. . .) denotes a Boltzmann ther- mal average. 
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The Czz, Cxx and Cx coefficients for Ho(OH)3 are plot- 
ted in Fig. [SI The inset shows a comparison of the 
terms in TLmf- In Ho(OH)3, the coefficient LxxCx^ (the 
fourth therm of TYmf) does not exceed 1.5% of tlie ef- 
fective transverse field, F = CxxCxCxx — ^^{Bx)- In 
Tb(0H)3, this ratio is even smaller, and we thus neglect 
it, further motivated by the fact that doing so decou- 
ples rriz from nix and make the problem simpler. The 
"^iriyay in Eq. [12] can be omitted, since due 



to symmetry my 



0. The interaction correction. 



t^xx Cx Cxx, to the effective transverse field, F, is of order 
of 3% of F and we retain it in our calculations. Thus, we 
finally write 



(13) 



where P = C^^ C^z'^m^ and F = Cxx Cx Cxx ~ \^{Bx)- 

Diagonalizing the Hamiltonian (fTSll allows us to eval- 
uate mz and nix = (c^), giving well known formulae^: 



nix 
ni. 



^^tanh(V P^TT^/ r), 
^^tanh(VP^TT^/T) ' 



and the phase boundary, 
TciBx) - 



TiBx) 



atanh( ^^^^j ) 



(14) 



(15) 



In Fig. [6l we show that Eq. p5|) yields a phase diagram 
that only insignificantly differs from the one obtained 
from the full diagonalization of TYmf in Eq. ([7]) shown 
in Fig. [21 in the case of IIo(OII)3, and, in the case of 
Tb(0H)3, the discrepancy is even smaller because the 
energy gap to the third crystal field state, 118 cm~^ ~ 
170 K, is very large compared to = 5.59 K. 

As alluded to above, in the case of Dy(0H)3, a descrip- 
tion in terms of an effective Ising Hamiltonian method 
does not work because of the admixing between states 
of the two lowest doublets induced by the local mean- 
field that is proportional to (J^) (see Appendix A). The 
dashed line in the last panel of Fig. [6| shows the incor- 
rect phase diagram obtained for Dy(0II)3 obtained us- 
ing an effective spin-1/2 Hamiltonian constructrcd from 
only the ground doublet. It turns out that a form of the 
method of Section HVl can still be used. However, instead 
of keeping only two levels in the interaction Hamiltonian, 
one needs to retain at least four states. In analogy with 
the procedure in Section [TVl we diagonalize the single ion 
Hamiltonian Tio of Eq. ^ which consist of the crystal 
field Hamiltonian and the transverse field term. Next, we 
write an effective interaction Hamiltonian using the four 
(or six) lowest cigcnstatcs of Hq. The resulting effective 
Hamiltonian is then used in the self-consistent Eqs. ([S]). 
For example, for Bx = 4.8 T, proceeding by keeping only 
the four lowest eigenstates of Ho to construct the effec- 
tive Hamiltonian, one finds a critical temperature that is 
only about 3% off compared to a calculation that keeps 
all 16 eigentates of Hq. This difference drops below 1% 
when keeping the 6 lowest cigcnstatcs of Tio . 




Figure 6: Comparison of the phase diagrams obtained with di- 
agonalization of the full manifold (black, solid lines) and with 
effective spin-1/2 Hamiltonian (red, dashed lines). Calcula- 
tion for Ho(OH)3 were performed using the CFP of Karmakar 
et alM^. 



Having explored the quantitative validity of the 
spin-1/2 TFIM description of Ho(OH)3 and Tb(0H)3 in 
nonzero Bx, we now turn to the problem of the first order 
PM to FM transition at large Bx and low temperature in 
Dy(0H)3, exposed in the numerical solution of the self- 
consistent equations comprised in Eq. Q (with ii = x, z). 



V. FIRST ORDER TRANSITION 

The first order transition in Dy(0H)3 takes its origin 
in the sizeable admixing among the four lowest levels in- 
duced by the the local mean-field that is proportional to 
(Jz). Under the right temperature and field conditions, 
two free-energy equivalent configurations can exist: an 
ordered state with some not infinitesimally small mag- 
netization, (Jz) > 0, and a state with zero magnetic 
moment, (Jz) = 0. To simplify the argument, we con- 
sider how this occurs at T = 0. At first, let us look at 
the situation when the longitudinal internal mean field 
induces an admixing of the ground state with the first 
excited state only (as in the TFIM). In such a case, there 
is only a quadratic dependence of the ground state en- 
ergy on the longitudinal mean field, B^^ , and, conse- 
quently, only one energy minimum is possible. Now, if 
there is an admixing of the ground state and at least 
three higher levels, the dependency of the ground state 
energy on B^^ is of fourth order and two energy minima 
are, in principle, possible. Thus, at a certain value of 
external parameters the system can acquire two energet- 
ically equivalent states, one with zero and the other with 
a non-zero magnetization. When passing through this 
point, either by varying the transverse field or the tem- 
perature, a first order phase transition characterized by a 
magnetization discontinuity occurs. To make this discus- 
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sion more formal, we now proceed with a construction of 
the Ginzburg-Landau theory for Dy(0H)3 for arbitrary 
Bx in the regime of and T values where the param- 
agnetic to ferromagnetic transition is second order. This 
allows us to determine the the tricritical transverse field 
value above which the transition becomes first order. 



A. Ginzburg-Landau Theory 

To locate the tricritical point for Dy(0H)3, we per- 
form a Landau expansion of the mean-field free energy, 
^mf{{Jx) ,{Jz))- Next, wc minimize J-'mf with respect 
to (Jx) , leaving (J^) as the only free parameter. The 
mean-field free energy can be written in the form 

^mf((J.),(J.)) = -T log Z{{Jx),{Jz)) 

- ^ (Cxx{Jxf +L,,{J,f^ ,(16) 

where Z {{Jx) , {Jz)) is the partition function. 

Just below the transition, in the part of the phase di- 
agram where the transition is second order, {J^) is a 
small parameter (i.e. has a small dimensionless numeri- 
cal value). We therefore make an expansion for {J^) as 
a function of (Jr), which we write it in the form: 

{Jx)^{Jx)o + 5{{Jz)). (17) 

(>^;)o is the value of {J.^) that extremizes Tmf when 
(Jz) = 0. 5{{Jz)) is a perturbatively small function of 
{Jz), which we henceforth simply denote 6, and which 
is our series expansion small parameter for {Jx)- Sub- 
stituting expression (fT7|) to Hmf of Eq. and setting 
Jcx = for the time being, we have 

= 'Hcf{Ji) — QI-IbBxJx 

+ CxxJx{{Jx)o+S)+C,,J,{Jz), (18) 

or 

n = no{Bx, {Jx)o) + ^xxJxS + c,,j, (J,) , (19) 

where for brevity, as in Eq. ([7]), the constant term has 
been dropped because, again, it does not affect the ex- 
pectation values needed for the calculation. 

The power series expansion of the partition function, 
and then of the free energy (fT6|) . can be calculated from 
the eigenvalues of Hamiltonian . Instead of applying 
standard quantum-mechanical perturbation methods to 
Eq. (fT9|) . we obtain the expansion of energy levels as a 
perturbative, 'seminumerical', solution to the character- 
istic polynomial equation 

det [Ho + CxxJxS + C,,J, {J,) - En] = 0. (20) 

We can easily implement this procedure by us- 
ing a computer algebra method (e.g. Maple™ or 



Mathematica"^^). To proceed, we substitute a formal 
power series expansion of the solution 

En = -f (Jz) {J.)'S + . . . , 

(21) 

to Eq. (|20p . containing all the terms of the form 

-Bn"''^' {Jz)"S^, where a 2/3 < 6, as wih be justified 
below Eq. ([M)) . To impose consistency of the resulting 
equation obtained from Eq. ([20]) and Eq. ([2T|) , up to sixth 
order of the expansion in {Jz) in), we need to equate to 
zero all the coefficient with the required order of {Jz) 
and (5, i.e. a + 2(3 < 6. This gives a system of equations 
that can be numerically solved for the coefficients E'^'"''' , 
where fc, / > 0. By E^'°^ we denote the eigenvalues of 
the Hamiltonian TCo{Bx, {Jx)q)- 

We use the perturbed energies, i?„, of Eq. (pij) to cal- 
culate the partition function 

Z(5,(Jz))=^e^^"/^ (22) 

n 

and substitute it in Eq. p^ . We Taylor expand the 
resulting expression to obtain the numerical values of the 
expansion coefficients in the form 

(23) 

The free energy J-'mf is a symmetric function of (Jz), so 
the expansion contains only even powers of (Jz)- We 
minimize J^mf in Eq. (j23[) with respect to S. To achieve 
this, we have to solve a high order polynomial equation 
dJ-MF/dS = 0. Again, we do it by substituting to the 
equation a formal power series solution 

6{{Jz))^D2{Jz)^ + Di{Jz)^ + ... (24) 

and then solve it for the values of the expansion param- 
eters Dn- Due to symmetry, only even powers of (Jz) 
are present and, from the definition of S, the constant 
(Jz) -independent term is equal to zero. From the form 
of the expansion in Eq. , we see that to finally obtain 
the free energy expansion in powers of (Jz), up to ??-th 
order, we need to consider only the terms (Jz)" where 
a + 2(3 < n. Finally, by substituting S from Eq. into 
Eq. (|23p . we obtain the power series expansion of the free 
energy in the form: 

.^MF = Co + C2 ( Jz)' + C4 ( Jz)' + Ce ( Jz)' . (25) 

In the second order transition region, the condition 
C2 = with C4 > parametrizes the phase boundary. 
The equations C2 = C4 = gives the condition for the 
location of the tricritical point. In the regime where C4 < 
0, the condition C2 = gives the supercooling limit. The 
first order phase transition boundary is located where 
the free energy has the same value at both local minima. 
Increasing the value of the control parameters, T and 
Bx, above the critical value, until the second (nontrivial) 
local minimum of J^viF vanishes, gives the superheating 
limit. 
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Figure 7: Tricritical behavior of Dy(OH)3. The continuous 
line marks the phase boundary. The dot indicates the location 
of the tricritical point. The upper and lower lines (dashed) 
are the hmits of the superheating and supercooling regimes, 
respectively. As an example, the inset shows the order pa- 
rameter (J,) vs Bx, at the temperature T —0.3 K. For this 
temperature, the phase transition occurs at Bj, ~ 4.98 T. 
The upper dashed extension of the solid line corresponds to 
the superheating limit. The dash-dotted line to the left of 
the tricritical point shows the transition at the supercooling 
limit. 



The location of the tricritical point is Tj^^ =0.75 K, 
BjCP =4.85 T. We show in Fig. [7] the first and the 
second order transition phase boundary; the tricritical 
point is marked with a dot. In the first order transition 
regime, the superheating and supercooling limits are also 
plotted. (Jz) ceases to be a small parameter for values 
of T and 'away' from the tricritical point. Thus, the 
two upper curves in the phase diagram of Fig. [7] are de- 
termined from a numerical search for both local minima 
of the exact mean-field free energy in Eq. p6)) without 
relying on a small (J^) and S{{Jz)) expansion. The su- 
percooling limit is calculated from the series expansion 
((25)) and determined by the condition C2 = 0. 

In the inset of Fig. [71 we show the average magnetic 
moment, (Jz), as a function of the transverse field, at 
the temperature of 0.3 K. The dots and the dashed lines 
mark the supercooling limit, first order phase boundary 
and the superheating limit, in order of increasing B^. 
The shape of the free energy at these three characteristic 
values of the magnetic field, Bx, at temperature of 0.3 K, 
is shown in Fig. [51 

In Fig. [51 wc plot the free energy as a function of (Jz), 
where (J^) is minimizing J-mf as a function of (Jj) at 
T = 0.3 K. Free energy at the phase transition (Bx ~ 
4.977 T) is plotted with a continuous line. The dashed 
and dot-dashed plots show free energy at the superheat- 
ing and supercooling limits, at Bx « 4.995 T and Bx ~ 
4.940 T, respectively. The free energy clearly shows the 
characteristic structure (e.g. barrier) of a system with a 
first order transition. It would be interesting to invcsti- 



Figure 8: Free energy vs average magnetic moment (Jz) at the 
temperature of 0.3 K. The free energy at the first order phase 
transition is plotted with continuous line. The upper (dashed) 
line shows the free energy at the boundary of the superheating 
limit. The lower (dot-dashed) line shows the free energy when 
the system passes through the supercooling limit. All the 
plots show the free energy displaced by a constant chosen 
such that their shape can be compared. 



gate whether the real Dy(0II)3 material exhibits such a 
Ba;— induced first order PM to FM transition at strong 
Bx ■ In the event that the transition is second order down 
to T = and Bx = B%, Dy(0H)3 would oflFer itself as 
another material to investigate transverse field induced 
quantum criticality (see 4'^ item in the list at the begin- 
ning of Section IV). However, a quantitative microscopic 
description at strong Bx would nevertheless require that 
the contribution of the lowest pairs of excited crystal field 
states be taken into account. 

One may be tempted to relate the existence of a first 
order transition in Dy(0H)3, on the basis of Eq. (j23p . 
with two expansion parameters {Jz) and 5, to the familiar 
problem where a free-energy function, J^{m,e), of two 
order parameters m and e. 



2 2 
— e — qem 
2 ^ 



displays a first order transition when g^/K > 6/2. How- 
ever, we have found that this analogy is not useful and 
the mcchamism for the first order transition is not triv- 
ially due to the presence of two expansion parameters, 
(Jz) and 6, in the expansion ([23]) . It is rather the com- 
plex specific details of the crystal field Hamiltonian for 
Dy(0H)3 that are responsible for the first order tran- 
sition. For example, at a qualitative level, a first order 
transition still occurs even if d{{Jz)), in Eq. ([T7]), is taken 
to be 0, for all values of (Jz)- 
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Figure 9; Average magnetic moment, (J,), vs transverse mag- 
netic field, Bx, ai T = 0, Jcx = 0, for different values of 
longitudinal magnetic field, B^- 



B. The effect of longitudinal magnetic field and 
exchange interaction on the existence of first order 
transition in Dy(OH)3 




Figure 10: Temperature corresponding to the TCP as a func- 
tion of nearest-neighbor exchange constant Jcx- The upper 
inset shows position of the TCP in the phase diagram plane 
for the same range of exchange constant as in the main plot. 
For an exchange Jcx > Jcx"*, with Jcx"* = 0.995, the tricritical 
point at T > ceases to exist. The lower inset shows average 
magnetic moment {J.) vs transverse field B^ at temperature 
T —0, for the shown values of exchange constant Jcx. 

Having found that the PM to FM transition may be 
first order in Dy(0H)3 at large (low T), it is of inter- 
est to investigate briefly two effects of physical relevance 
on the predicted first order transition. Firstly, since the 
transition is first order from < T < rTCP ^ may ask 
what is the critical longitudinal &c\d, B^, required to push 
the tricritical point from finite temperature down to zero 



netic field, Bz, applied along the longitudinal z direction 
destroys the first order transition, giving rise to an end 
critical point. We plot in Fig. [9] the magnetization, (Jz), 
as a function of B^ for different values of i?^ at T = 0. 
We see that a critical value of Bz is reached between 1 
T and 2 T, where the first order transition disappears, 
giving rise to an end critical point at T = 0. Hence, as- 
suming that the low-temperature B^-driven PM to FM 
transition is indeed first order in Dy(0H)3, the results of 
Fig. ^ indicate that the critical longitudinal field for 
a quantum critical end point is easily accessible, using 
a so-called vector magnet (i.e. with tunable horizontal, 
Bx, and vertical, Bz, magnetic fields}^. 

It was discussed in Section llIll (Fig. [4|) that the (yet un- 
determined) nearest-neighbor exchange interaction, Jox, 
affects the zero B^ critical temperature, Tc, and the zero 
temperature critical transverse field, B^ . It is also of in- 
terest to explore what is the role of Jox on the location 
(temperature and transverse field) of the tricritical point 
in Dy(0H)3. 

We plot in Fig. [10] the temperature corresponding to 
the tricritical point (TCP) as a function of antiferromag- 
netic exchange and, in the upper inset, the location of the 
TCP on the phase diagram is presented. The location of 
the TCP was calculated using the semi-analytical expan- 
sion described in Section IVAI Wc found that the system 
ceases to exhibit a first order transition at nonzero tem- 
perature when the value of nearest neighbor exchange 
constant, Jox, exceeds Jo"'^ = 0.995. At Bx=0, the crit- 
ical temperature calculated with the value of exchange 
constant Jox=0.995 is 4.09 K. In the lower inset of Fig. 
[TOlwc plot the average magnetic moment, (Jz), as a func- 
tion of Bx at zero temperature, for different values of Jox. 
The top inset shows a parametric plot of the position of 
the TCP in the (T, B^) plane as Jcx is varied. 



VI. CONCLUSION 

We have presented a simple mean-field theory aimed 
at motivating an experimental study of transverse-field- 
induced phase transitions in the insulating rare-earth 
Ising RE(0H)3 (RE=Dy, Ho) uniaxial dipolar ferromag- 
netic materials. 

In setting out to perform the above calculations, we 
were mostly motivated in identifying a new class of ma- 
terials as analogous as possible to LiHo2;Yi_j^F4, where 
interesting phenomena, both in zero and nonzero applied 
transverse field B^, have been observed. In particular, 
we were interested in finding compounds where a sys- 
tematic comparison between a non-Kramers (e.g. Ho"^"*") 
and a Kramers (e.g. Dy'^+) variant could be investigated. 
From our study, we are led to suggest that an experi- 
mental study of the Dya;Yi_2:(0H)3 and Ho2;Yi_2;(OH)3 
materials could bring new pieces of information on the 
physics that may be at play in LiHoxYi_2;F4 and to as- 
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certain if that physics is unique to LiHoa;Yi_a;F4 or if it 
also arises in other dihited dipolar Ising ferromagnets. 

Depending on the details of the Hamiltonian charac- 
terizing Dy(0H)3, it may be that a first order transition 
occurs at low temperature (large B^), due to the admix- 
ing between the ground doublet and the low-lying crystal 
field states that is induced by the spin-spin interactions. 
For the same reason, we find that Dy(0H)3 is not well de- 
scribed by an effective microscopic transverse field Ising 
model (TFIM). On the other hand, Ho(OH)3 appears to 
be very well characterized by a TFIM and, therefore, con- 
stitutes a highly analogous variant of LiIIoF4. Tb(0II)3 
is also very well described by a TFIM. Unfortunately, in 
that case, the critical B^, B"^, appears prohibitively large 
to be accessed via in-house commercial magnets. 

We hope that our work will stimulate future system- 
atic experimental investigations of these materials and, 
possibly, help shed some light on the rather interesting 
problems that pertain to the fundamental nature of clas- 
sical and quantum critical phenomena in disordered dipo- 
lar systems and which have been raised by nearly twenty 
years of study of LiIl0a;Yi_xF4. 



Acknowledgments 

We thank P. Chakraborty, R. Cone, L. Corruccini, J.- 
Y. Fortin, S. Girvin, R. Higashinaka, R. Hill, Y. Maeno, 
B. Malkin, H. Molavian, V. Oganesyan, T. Rosenbaum, 
D. Silevitch, A. Tabei, K.-M. Tam, W. Wolf, and T. Ya- 
vors'kii for useful discussions. This work was supported 
by the NSERC of Canada, the Canada Research Chair 
Program (Tier I, M.G), the Canadian Institute for Ad- 
vanced Research, the Canada Foundation for Innovation 
and the Ontario Innovation Trust. M.G. thanks the Uni- 
versity of Canterbury (UC) for funding and the hospital- 
ity of the Department of Physics and Astronomy at UC 
where part of this work was done. 



Appendix A: PERTURBATIVE CALCULATION 
OF THE PHASE DIAGRAM IN DY(OH)3 

To investigate the role of the Jz matrix elements be- 
tween the two lowest states and the first excited levels 
on the magnetic behavior of Dy (011)3, we calculate the 
critical temperature for a second order transition using 
second order perturbation theory. This method is exact 
in the second order phase transition regime, where B^ is 
less than the tricritical field value, BJ^p (BJCP ^ 4 35 
when using the CFP of Scott et al.). 

For a given value of the transverse field, B^, and the 
corresponding value of average magnetization in trans- 
verse direction, {Jx), we consider CzzJz {Jz) term as a 
perturbation to the reference mean-field Hamiltonian, 



'Ho — Hcf{Ji) — gf^sBxJx + CxxJx {Jx) , 



(Al) 



describing the PM phase at a temperature T > Tc{Bx) as 
in Eq. [7]Here too, we have dropped the constant terms. 
The eigenvalues, Ep, and eigenstates, \p), of the per- 
turbed Hamiltonian, 



H =Ho + CzzJz {Jz 



(A2) 



are written in terms of eigenvalues, Ep^\ and eigenstates, 
Ipf^")), of the unperturbed Hamiltonian, Hq, of Eq. (|A1[) . 

Ep = + {Jz) + {Jz)' E^^\ (A3) 

b) = b(°)) + (J.)5]c(^i|fcW). (A4) 

The coefficients of the pcrturbative expansion are given 
by 



J^p — '-'ZZ'Jpp, 



and 



^pfe 



E 



^'zz\Jkp\'^ 



E 



(0) 



E 



(0) 



'kp 



E; 



(0) 



E 



{0)- 



(A5) 



(A6) 



(A7) 



where J|p = (fc(°) | J. are the matrix elements of 

the Jz operator in the basis of eigenvectors of the unper- 
turbed Hamiltonian, Tio- The applied magnetic field, Bx^ 
lifts the degeneracy of the Kramers doublets, thus we can 
use the nondegenerate perturbation method. The diag- 
onal elements of the Jz operator vanish, hence, the first 
order correction to energy vanish, iJp^' = Jpp = 0. 
We calculate the thermal average of the Jz operator, 



{Jz) 



V 



{p\ Jz \p) e 



/z, 



(A8) 



using the perturbed eigenstates, \p), and eigenvalues, Ep, 
and where Z = ^pe~^^^'^ . Keeping only terms up to 
third order in {Jz) in the expansion of Eq. (jA8|) . we find 



z 



4")(l + ifp(J,)') (A9) 



and 



{p\ Jz \p) = 2Czz {Jz) 



\Jpk? 



(0) 



k^p 



where, for convenience, we write n 



E, 



(0) ■ 



zc 



0) 



(0) 
p 



(AlO) 



and K„ 



Thus, we can write 
(J.) = 2(J,)^nf (l + Kp{Jz)')Y. 



Czz\J. 



pk\ 



p(0) p(0) ' 



k^p ^k 

(All) 
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and finally, we get 

(J. 



2 _ ^l^pM^p'^P £;(»)-£;<«) 



2V ^"'^g^'' 



(A12) 



Putting (J^) = in Eq. ((IT2l) wc obtain the condition 
for the critical temperature Tq. 



2E 



(0) R.(0) 



- E 



(A13) 



fe 



In solving Eq. (jA13[) for Tc, we have to self-consistently 
update the value of {Jx) in order to diagonalize Tio in 
Eq. dJT]) and to find . Solving Eq. ((Xl3l) with only 
the 4 lowest energy eigenstates (after diagonalizing the 
full transverse field Hamiltonian of Eq. ([5])), yields a 
phase diagram that is in good agreement for < Sj^^ 
with the phase boundary obtained with all crystal field 
eigenstates (or equivalently from Eq. ([H])). 

Estimating the values of the elements of the sum in 
Eq. (jA13[) . one can see that the matrix elements of the Jz 
operator, mixing the two lowest states with the excited 



states, may bring a substantial correction to the value 
of the critical temperature obtained when only the two 
lowest eigenstates are considered. In the low temperature 
regime, one could omit the matrix elements between the 
states of the excited doublet, but we have to keep the 
matrix elements between the states of the ground dou- 
blet and first excited doublet. The contribution from the 
further exited states is quite small because of the increas- 
ing value of the energy gap present in the denominator 
in Eq. (^13)) . 

At r=0, Eq. (|A10[) leads to the equation for the critical 
transverse field, B'l: 



Lzz\JUB%)^^ 



Efm 



(A14) 



Again, we see that the matrix elements of operator, 
admixing the ground state with excited levels, have to 
be considered. Note that since Eq. (|A14[) pertains to the 
case of zero temperature, this equation is only valid in 
a regime where the transition is second order (i.e. when 

^cx ^ •■'cx )■ 
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